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In this paper we have considered a Gaussian polymer chain of length L as an intrinsic object 
enclosed on a surface embedded in the Euclidean space. When the surfaces are the sphere S 13-1 in 
D dimensions and the cylinder, the cone and the curved torus in R 3 , we have calculated analyti- 
cally and numerically (using the diffusion equation and the path-integral approach) the probability 
distribution function G(R|R'; L) of the end-to-end vector R — R' and the mean-square end-to-end 
distance ((R — R') 2 ) of the polymer chain. Our findings are: the curvature of the surfaces induces a 
geometrical localization area; at short scales (L — » 0) the polymer is locally flat and the mean-square 
end-to-end distance is just the Gaussian value LI (I — Kuhn length), independent of the metric 
properties of the surface; at large scales (L — > oo), ((R — R') 2 ) tends to a constant value in the 
sphere case, it is linear in L for the cylinder and reaches different constant values (as a function of 
the geometry of the surface) for the curved torus. In the case of the cone, contraction of the chain 
is induced at all length scales by the presence of the vertex, as a function of the opening angle 2a 
and the end position R' of the chain. Explicit crossover formulas are derived for G(R|R';L) and 
((R — R') 2 )- 



I. INTRODUCTION 



Statistics and topological characteristics of polymers 
embedded in curved interfaces are of significance in bi- 
ological problems related to the way polymers wrap 
around spherical and rodlike macromolecules and vesi- 
cles, and the role played by geometry and topology in 
the dynamics of a biomolecule captured on the mem- 
brane of a cell These types of issues arise in other 
area of physics as well: the quantum mechanics of free 
particles in curved and multiply connected spaces (e.g. 
see Ref. 0), the physics of vortex lines in Type II super- 
conductors H and the rotational Brownian motion. 

In realistic situations, excluded-volume interaction be- 
tween polymer segments and fluctuations of curved in- 
terfaces need to be accounted for. Under these circum- 
stances, only approximate analytical results can be gener- 
ally obtained via perturbation analysis, variational tech- 
niques or renormalization group theory. Before embark- 
ing on such an effort, it is necessary to solve the problem 
of a free polymer chain, characterized only by its connec- 
tivity, on a curved, fixed surface. 

In the present paper we derive, using path-integral 
and spectral representation, the probability distribu- 
tion of the end-to-end vector G(R|R'; L) and the 
size — expressed by the mean-square end-to-end distance 
((R - R') 2 ) in the ambient space — of an ideal (non- 
interacting) Gaussian polymer lying on a curved man- 
ifold. 



The curved surfaces studied in this paper are the S 
sphere in D dimensions (Fig. |l|A, for D — 3); the cylin- 
der (Fig. [pB), the circular cone (Fig. ^C) and the torus 
(Fig. |I]d) in R 3 . The geometrical parameters are shown 
in Fig. |l|. The polymer is an intrinsic D — 1-dimensional 
object in the spherical surface, while it is 2-dimensional 
in the other surfaces. 

Although, due to the equivalence between the random- 
walk problem and the statistics of an ideal Gaussian chain 
Q , many results are already available from the theory of 
Brownian motion, the specific problems we address here 
have not been answered yet. 

In our calculations we have used the path-integral rep- 
resentation of the probability distribution function 
G(R|R';L) of a Gaussian polymer chain and the diffu- 
sion equation obeyed by this probability. As mentioned, 
no excluded volume or other types of interactions are ac- 
counted for in this paper. 



II. FORMULATION 



Let us consider a D-dimensional Gaussian polymer 
chain with N links of Kuhn-lcngth I and total length 
L — Nl, with the first bead located at R' and the last 
one at R. The end-to-end vector R — R' is distributed 
according to the following probability distribution func- 
tion, given as an explicit path integral || 
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where D is the intrinsic dimension of the polymer, 5s 
is the change in the arc-length s along the chain and 
Yj is the position vector of the j-th bead of the polymer 
chain. Note that the distance on the surface is calculated 



as {vj 



= 9uv(Uj 
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-i), where u, v 



are the intrinsic coordinates and the metric tensor g uv 
is discretized in a symmetric form or using the midpoint 
rule (e.g. see Ref. [ HQ]). In R D , the expression (2.1) 
gives the Gaussian probability distribution: 
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Here M and are obtained from the eigenvalue equa- 
tion 
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with fi some generalized index. In general, the probabil- 
ity distribution G(R|R'; L) is not translational invariant. 
The mean-square end-to-end distance can be calculated 
from 
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It is to be noted that the calculated mean-square distance 
represents the distance between the ends of the polymer 
measured by an outside observer, located in the ambi- 
ent space, which is different from the geodesic separation 
measured along the surface. 



The probability G(R|R'; L) obeys the diffusion equation 
^-±a)g ( R|R';L)=O i (2.3) 



subject to the initial condition 

lim G(R|R';L) = 5 (D) (R-R'). 



(2.4) 



When the polymer lies in a curved manifold, one needs 
to employ the Laplace-Beltrami operator given by (e.g. 
see Ref. 
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III. CALCULATIONS AND RESULTS 



Due to the kn own equivalence [j4||8j between the dif- 
fusion equation fl2.3| ) and the Schrodinger equation for 
a free particle, we can readily obtain the polymer prob- 
ability distribution G(R, R';L) for the sphere and the 
cylinder cases, from the corresponding probability am- 
plitudes for a rigid rotor in D dimensions |^] and for a 
2-D free particle moving on the surface of a cylinder, 
respectively. 



A. The Sphere S D1 in 



where is the inverse of the metric tensor g^ of the 
surface and g = det(<?jj). We remark that gij is the in- 
duced metric from the ambient embedding space (R D or 



One can obtain the solution of Eq. (|2.3| ) satisfying the 
initial condition (2.4) by applying a Laplace transform 
with respect to L that gives the Green's function equa- 
tion 



E-^A\G(R\K';E) = 5( D \R-R'), (2.6) 



which can be solved by an eigenfunction expansion of G 
followed by an inverse Laplace transform that eventually 
yields the heat kernel: 



Let a be the radius of the D — 1 dimensional sphere 
centered at the origin of the coordinate system. Then 
R(s) is a vector of magnitude a that measures the posi- 
tion of the bead at s about the center of the sphere. The 
intrinsic dimension of the polymer is D — 1. In spherical 
coordinates the Laplace-Beltrami operator ( [2.5] ) becomes 
§ §4.2.3] 



A = \c\ 



(3.1) 



where C 2 is the Z3-dimcnsional Legendre operator, satis- 
fying the eigenvalue equation 



£ 2 Sy (fi) = -J(J + D — 2)S?(fi). 



(3.2) 



G(R|R';i) = E^ I (R)0*(R')exp(-^ 1 i). (2.7) The Sj n functions are the hyperspherical harmonics or- 
M thonormal and complete on the unit S " 1 sphere JlO| | : 
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where M — ( 2J + g ^^) (■^+- p ~ 3 ) ! an( \ fi [ s fch e solid angle 

in D di mens ions. Using Eq. ( |3.2|) and the eigenvalue ex- 
pansion Q, where </y -> S} 1 and -> J(J + D - 2), 
we readily find the probability distribution function of a 
polymer chain (or equivalently, of a random walk) that 
starts at R' and ends at R in TV = L/l steps, on the 
S ^ 1 sphere: 
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The probability found above is normalized properly to 
one (J d D R G(R|R; L) = l). In K 3 we recover the known 
lOl rotational diffusion result: 
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Because G(R|R'; L) is translationally invariant, the end- 
to-end distance can be calculated as a conditional ex- 
pectation value by dropping the integration over R' in 
Eq. (|]t]). To evaluate the remaining integral, we use the 
addition theorem for hyper-spherical harmonics p~C| ] 

f; sy(n)sr(nO = ^ (2J ^ 2 ' 2) 



rn — 1 



x C 



(D-2)/2 
J 



(cos*rr<), (3.6) 



9 D / 2 

and Cj are the Gegenbauer poly- 



where Sd 

normals of argument cos?/>rr'. Writing (R — R') 2 
2a 2 [1 — cos^rr/] and inserting (3.4) in (|2.9D, we get 
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Choosing R' along one of the coordinate axes and apply- 
ing the recurrence relation (lj, §8.933] 



cos Of (cos 6) = + !) CT+i(«»ff) 

+ (2i/ + /-l)Cf_ 1 (coBfl)], (3.8) 

we notice that the orthogonality of the Gegenbauer poly- 
nomials (with Cq(cos8) = 1) imp lies that all terms ex- 
cept J — 1 are zero. From ([3/7]) we finally obtain the 
mean-square end-to-end vector for a polymer living on 
the surface of a D — 1-dimensional sphere as: 
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Therefore, a characteristic area 2a 2 called geometrical lo- 
calization area emerges in the description of ( (R — R') 2 ) > 
due to the geometrical restriction of space. 

In contrast to the problem of a stiff polymer chain 
without any confinement, where the chain expands due 
to chain stiffness, chain contracts here by geometrical 
confinement. If the area of the confining sphere is small 
in comparison with LI, then the polymer wraps many 
times around the sphere so that ((R — R') 2 ) approaches 
the geometrical localization area. On the other hand if 
a 2 3> LI, then the polymer does not sense the curva- 
ture of the sphere so that ((R — R') 2 ) — LI. The mean- 
square end-to-end distance (3.9) normalized at a 2 is plot- 
ted in Fig. | as a function of the dimensionless variable 
L\. It is also to be noted that ((R — R') 2 ) is independent 
of D. 



B. The Cylinder in 



Consider a cylinder of radius a and of infinite length, 
with the induced metric from R 3 (in cylindrical coordi- 
nates, with Oz along the axis of the cylinder): = a 2 ; 
94>z = 0; g zz = 1; ^fg = a. The position vector R(s) of 
any bead s has the components (a cos </>, a sin </>, z). The 
polymer dimension is D = 2. 

The Green's function equation in the Laplace space 



(2.6) can be solved directly by expanding the G function 
in the complete set of orthogonal functions on the circle 
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and Fourier transforming with respect to z - 
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After applying an inverse Laplace transform and eval- 
uating the integral over fe, the spectral expansion of 
G(R|R';i) is found to be 
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An equivalent form of the expression above can be ob- 
tained by using the relation 
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and the Poisson formula 
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Applying these transformations in Eq. ( |3.11 ), we arrive 
at the the winding number expansion 
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where n (a topological term) is the winding number that 
counts how many times the chain winds around the cylin- 
der. A term in the sum represents now the probability 
that a given chain starting at R and ending at R' will 
wind around the cylinder n times (in L/l steps). This 
expansion can be used to investigate the statistics of the 
winding number (e.g. the mean-square winding number 
(N 2 )). 

To calculate the mean-square end-to-end vector, we 

2q 2 [1 — cos(0 — (j)')] and 
in (2.9). Evaluating the 
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integrals, we get 
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As expected, the polymer remains almost Gaussian (a 
prefactor of 1/2 instead of 1) along the Oz direction. 
The influence of the geometry appears through the ra- 
dius a, related to the mean curvature of the cylinder. 
Once again, the polymer becomes geometrically localized. 
The mean-square end-to-end distance (normalized at o 2 ) 
is represented in Fig. |as a function of the dimensionless 
variable ^4. 



C. The Cone in 



with the opening angle 2a. In cylindrical coordinates, 
the position vector R(s) of any point s along the chain 
is R(s) = (p{s) cos (j)(s),p(s) sin (f>(s),p(s) cot a) and the 
induced metric is: g pp = 1 + cot 2 a; g p( j, = 0; g^ = p 2 ; 
■Jg = pVl + cot 2 a. As in the previous problem, the 
dimension of the polymer is D = 2. 

To compute the probability distribution G(R|R'; L) of 
the end-to-end vector we start with the Green's function 



equation (2.6) obeyed by the Laplace transform G: 
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We solve this equation in a standard manner [[L3|. First 
we eliminate the azimuthal dependence of G by inserting 
the expansion 
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in (BTq), which gives 
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where we introduced the new variable x — kp, with 
k 2 = , A \ and v 2 — .'"> . The solutions of the ho- 

/ sm^ a sin' a 

mogeneous equation are the modified Bessel functions 
I|„|(a;) and K|„|(a;). Then, the regularity of g v {x\x';E) 
at x = 0, x — > oo, the continuity at x = x' and the jump 
in the first derivative impose the solution 
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Here p< and p > denote the smaller and the larger, re- 
spectively, between the variables p and p'; we have also 
used the fact that the Wronskian W of the modified 
Bessel functions is W[l v , K J (x ) = —-. Inserting the so- 
lution found back in Eq. ( |3 -17 ) and performing the in- 
verse Laplace transform using the formula |l4| , §5.16.56] 

C- 1 {K y [(v^+ Vb)VE] l v [{Va- Vb)VE]\ = (3.20) 
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Now let us consider the surface of cone with the ver- 
tex at the origin O, centered about the Oz axis and 



we arrive at the exact spectral expansion of G(R|R'; L) 
for a polymer on a cone: 
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which is properly normalized to one (JdRG(R|R'; L) — 
1). For a = I we recover the expression in polar coordi- 
nates of the propagator of a Gaussian polymer chain in 
the plane g||. 

As the total partition sum gives the area of the cone, 
we calculate the end-to-end distance by fixing R' and 
integrating only over R. Using cylindrical coordinates, 
with (R-R) 2 = (p 2 + p' 2 )/sm 2 a + 2pp'[l-cos((j)-(j)')} 
and applying various formulae Jl2[ involving definite in- 
tegrals of Bessel functions, we obtain 
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where t = p' /(Ll) 1 / 2 and 1F1 is the confluent hypergeo- 
metric function. 

Rather surprisingly, the mean-square distance depends 
on the position of the chain through the parameter t. To 
interpret this finding, we calculate the limiting values of 
((R — R') 2 ) for t — ■> 00 (polymer far on the conical sur- 
face or ct — ► 0), t — ► (polymer with one end attached 
to the cone vertex) and a = ^ (polymer on a 2—D plane 
surface). One finds that 

((R - R') 2 ) = Ll for {t = 0; t -> 00; a = -}, (3.23) 

which is just the Gaussian result in all cases. Physi- 
cally, when the polymer chain is either fixed at the origin 
(t = 0) or far away from it (t — > 00), there are no other 
preferred points on the cone and ((R— R') 2 ) has the 
planar value (we should also remark that a cone is ge- 
ometrically flat). In the crossover region, the singular 
character of the cone vertex becomes dominant and in- 
duces the contraction of the chain (possibly by winding 
about the origin). This behavior is illustrated in Fig. 
and Fig. ||, where the mean-square distance ((R — R')/ 
normalized at Ll is plotted as a function of the t param- 
eter and the angle a. 



D. The Torus in 



Consider a torus embedded in R 3 , with the circular 
cross-section of radius a. The axial circle containing the 
centers of the circular sections has the radius b > a (see 



Fig. ||). The torus is centered about the Oz axis and 
symmetric with respect to the xOy plane. We introduce 
the toroidal coordinates {r?, #,</>}, [|l5[ §2.13] where cf> is 
the usual azimuthal angle — in cylindrical coordinates — 
about Oz: 
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where < r\ < 00, < 6 < 2tt and < <f> < 2ir. The 
surface of the torus corresponds to a fixed value rjo of the 
coordinate r\. Then we have the following geometrical 
relations: 



b = c coth 770 ; a = — 



sinh 7/0 
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The induced metric tensor is g$Q 
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i r-7. We recall that the length interval on the 

(cosh r/o— cos v) z ° 

surface is ds 2 — ggg d9 2 + d<f> 2 and the volume element 
is d 2 r = y/gd9d(j). 

The calculation of the propagator of a Gaussian chain 
on a torus in R 3 (which is equivalent to solving the prob- 
lem of a free quantum particle moving on the same sur- 
face) cannot be done in closed form, as involves find- 
ing the eigenvalues of a Hill-type equation. Still, using 
the path-integral formalism and making certain "sloppy" 
approximations, we were able to derive an approximate 
form of G(R|R'; L) that behaves physically correct when 
computing moments of the distribution of R — R', which 
are the experimentally relevant quantities we are inter- 
ested in. This approach will be justified a posteriori by 
showing that one recovers the usual flat surface case and 
the expected asymptotic regimes obtained for spheres 
and cylinders. 



We start be recalling the path- representation (2.1) 
written in a compact form as 
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where the integration measure is given explicitly by 
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and Sj is the short-length polymer action 
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In the previous expressions, Ssj 



s i-i 



(with 



j = 1,N, so = 0, sjv = L) and Yj is the position vec- 
tor about the origin of the j-th element of the discretized 
chain. 

In order to evaluate the path-integral, due to the men- 
tioned equivalence between the statistical mechanics of 
polymers and quantum mechanics, we will apply the 
method used to calculate the propagator of a free particle 
moving on a circle and in the presence of a ring-shaped 
defect §,UlF 



We start with the action element Sj . Using the length 
interval on the torus, we write Sj in a symmetrized form 
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Observing that the integrals over the angular coordinates 
are similar to the path integral for a free particle moving 
on a circle (e.g. see Ref. [ PQ]), we have the formula 
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where {6j = 9(sj),(j)j = 4>{sj)} are the toroidal coor- 
dinates of the element located at Sj along the chain, 
A6j = 6j — 0,j_i and Atfij = <fij — <j)j-\. 

Next we apply a local rescaling which reads 

<jj = Ssj (cosh 770 — cos 9j ■) (cosh 770 — cos0 3 _i) (3.30) 

and which satisfies the global scaling relation: 

TV 

(7 = a j = i(cosh?7o — cos 9) (cosh 770 — cos 9'). (3.31) 
i=i 

This constraint is nontrivial, as it must be compatible 
with the local scaling (3.3C). That it gives the correct 
answer when applied in the path integral and when one 
uses an anisometric discretization (unequal length inter- 
vals Ssj), was discussed by Inomata |r7|| . Basically, it 
amounts to the prescription ]l8[ ] of ignoring terms of 
0(Ss^ +e ) in the path- integral, if e > 0. Then, the dis- 
cretized functional measure becomes: 

T)[r] — (cosh?7o — cos 9) (cosh 770 — cos 9') 

JV ^ JV-1 

x J[ — Yl d9j dfa c 2 sinh 770. (3.32) 
3=1° 3 3=1 

We apply now the length rescaling in the action Sj given 
by ( |3.29 ) directly, by dropping boldly any terms coming 
from a rigorous expansion of A9j and A<j)j about the new 
length variable Uj. As mentioned, this will be justified 
by the final results. Combining then with the integration 
measure, we get the approximate discretized propagator 
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Here a is an arbitrary positive constant and ipj = ipi^j) 
is the discretized angular variable on the circle, with ip' 
and i\) denoting the initial and the final positions. Equiva- 
lently, recognizing that the angular path integral contains 
trajectories that wind a different number of times around 
the origin, one can use the covering space mapping 



Returning to Eq. ( p.33| ) of the approximate discretized 
propa gator Gat, integrating over the angles usin g for mula 
( |3.34 ), and replacing a with its global value ( |3.31 ), we 
finally obtain — up to a normalization factor — the approx- 
imate distribution function of the end-to-end distance for 
a polymer chain on a curved torus in R 3 , expressed as a 
winding number expansion: 
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Another compact and illuminat ing f orm is obta ined by 
employing the transformations ( 3.12| ) and (3.13), which 
yield the angular momentum expansion 
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where we recall that cr(9, 9') — L(cosh?7o — cos#) (cosh 770 
— cos 9'), c 2 = b 2 — a 2 and cosh 770 = -. 63 is the theta 
function defined as |l^] 
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(3.37) 



G 



First, we stress that this is not the exact Green func- 
tion of the heat equation on the curved torus. Still, as 
necessary, the expression is symmetric in R, R', transla- 
tionally invariant with respect to the azimuthal angle cj>, 
and in the limit 770 — * 00 = a — ► reduces properly to 
the probability distribution of the end-to-end vec tor f or 
a polymer on a circle (see the angular part in Eq. (3.14)). 
Also, it obeys the initial condition 



G(R|R';0) = 5 (2) (R-R), 



(3.38) 



where N(x,t) is given by 
N(x,t) = 

x9 3 



<® dd ' 7 7^7 

(t-COS0)(f-COS0') 



' ix(t - cos0)(f - cos 6»') 



4vr(i 2 - 1) 



(3.44) 



The partition function Z{x 1 1) of the polymer on the torus 



and its trace has, at least numerically, the correct limit 
El when H -> 0: 



Z(x-+0,t) = JJd 2 RG(R,L\R,0) 



(3.39) 



, t 7T 3 t (t - Vt 2 ~ 1) , o, 

X 60 y/t 2 - 1 



Z(x,t) = a 2 (i 2 - l)*N(x,t). 



(3.45) 



This reduces properly to the area Air 2 ab of the torus when 
z> 1, t > 1 (6(z,t) ~ 1): 



Z(.r. i ) - 4tt 2 ab \jl - ^ ~ 4ir 2 ab (3.46) 



where we introduced the dimensionless variables: 

b LI 
t = cosh 770 = — ; x = . 



(3.40) 



A closed formula for ((R — R') 2 ) cannot be readily ob- 
tained and a numerical calculation is required. To begin 
with, we write the distance (in R 3 ) between the ends of 
the polymer chain: 



(R - R ) 



l\2 



2c 2 



(cosh 770 — cos 8) (cosh 770 — cos6>') 



— [cosh 770 



sinh 770 cos( 



?')]■ (3.41) 



Using the Green function (3.36) in the formula for the 
mean-square end-to-end distance (^^) with the proper 
volume element on the surface of the torus, evaluating 
the integrals over the azimuthal angles <f> and <j)' by ap- 
plying the definition (3.37) of the theta function and the 
formula: 



2jt 



d^d^ cos(0 - <t>') e"™(^'> = 27r 2 (5_i m + S+i m ), 

(3.42) 



and expressing all quant ities in terms of the dimension- 
less variables from Eq. (3.40), we eventually obtain the 
mean-square end-to-end distance as 



<(R-R') 2 ) _ 2(t 2 -l) 



t 2 + (l-t 2 )e 



N(x,t) J J 

E(t-OB8)(t-OB 0') 



dode' 



cos(6> - 6') 



x8 3 



(t - cos 9) 2 (t- cos 9') 2 
9' ix(t - cos 9)(t- cos 9') 



(3.43) 



4vr(t 2 - 1) 



Physically, for different values of x and t, three regimes 
are expected: 



a) LI < a 2 < b 2 

b) a 2 < LI < b 2 

c) a 2 < b 2 < LI 



t 2 



<x<l, 



— < K x, 
1 « ^ < x, 



(3.47a) 
(3.47b) 
(3.47c) 



where in the first case we should recover the solution for 
a polymer in the plane, in the second we must obtain 
the result for a polymer on a cylinder and in the third, 
the mean square end-to-end distance should reach a con- 
stant value (corresponding to a large winding number 
about both the Oz axis and the axial circle of radius b). 

These regimes are manifest in Fig. [|, where the re- 
sults for the torus are plotted with thin lines, for t — 
{2; 4; 6; 8; 10}. 



Initially (the first regime in Eq. (3.47)), the polymer is 



too small to explore the geometry of the surface, and we 
recover (numerically) the case of a polymer in a plane. 
This is also the behavior found for the other surfaces. 

As the length of the chain increases, it starts wind- 
ing along the circle of radius a and then along the axial 
circle, about the Oz axis. When t increases, the behav- 
ior is similar — for a certain range of x values — to that 
of a chain on a cylinder: the mean-square end-to-end 
distance is linear in LI/ 2. We can check a nalyt ically 
that the approximate propagator from Eq. (3.36) pro- 
duces the correct behavior. Because a 2 <C LI <C b 2 , the 
term n = dominates in the theta function in Eq. (3.43) 
(63(2, r) « 1) and the exponential coefficient of 1 — t 2 
can be expanded, so one obtains explicitly 



<(R-R') 2 ) ^2a 2 



y (a 2 « LI « b 2 ), (3.48) 
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which is just the li mit fo r the polymer on the cylinder of 
radius a from Eq. ( 3.15 ). 



Eventually, at any given t and large enough x, 
((R — R') 2 ) departs from the behavior of a polymer on 
a cylinder and approaches a constant value, which is a 
function only of the geometrical parameters a and b. This 
limiting value can be explicitly calculated. Now we have 
a 2 < b 2 -C LI, thus once again 63(2, r) » 1 but the ex- 
ponential coefficient of 1 — t 2 in the formula for distance 
( 3.43 ) is zero, which gives the limit: 

((R - R') 2 ) ~ 2 (a 2 + b 2 ) (a 2 <6 2 <L0- (3.49) 



If either a or b radii becomes 0, one recovers the asymp- 
totic limits of the end-to-end distance for a polymer on a 
circle or on a sphere, respectively, as found in Eq. ( |3.9| ). 

Although the asymptotic behavior is correctly recov- 
ered, we expect our results in the crossover region of in- 
termediate values of i and x to be only approximate. 



end-to-end distance of a Gaussian chain on the geomet- 
rical parameters of the confining surface. The for mulas 
for spheres and cylinders are simple (Eqs. (3.9, 3.15 )), de- 
scribing the crossover between free polymer regime and 
confined polymer regime. For the case of a cone, the fi- 
nal expression is more complicated, but still exact. For 
a torus, ther e are three asymptotic regimes as outlined 
by Eq. (3.47), and the crossover behavior is obtained nu- 
merically from Eq. ( 3.43] ). 

A number of open problems come into attention: the 
calculation of the mean-square geodesic end-to-end dis- 
tance (as measured on the surface and not in the em- 
bedding space), the statistics of the winding numbers 
(mainly for the polymer on the torus case), the influ- 
ence of the topological defects existing on the surface, 
the presence of excluded volume and other types of po- 
tential interactions. We hope to address some of these 
issues in the near future. 



IV. CONCLUSIONS 
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Due to its relevance to studies of the behavior of poly- 
mers confined at interfaces, we have considered here the 
problem of a linear, Gaussian polymer chain embedded in 
the following surfaces: the S ^ 1 sphere in D dimensions, 
the cylinder, the cone and the curved torus in R 3 . 

We obtained closed formulas for the probability dis- 
tribution function G(R|R';L) of the end-to-end vector 
of the chain for the sphere, the cylinder and the cone 
and an approximate propagator for the torus. We calcu- 
lated, analytically in the case of the sphere, the cylinder 
and the cone, and numerically for the torus, the mean- 
square end-to-end distance (in the embedding space) of 
the chain. As such, at least in the limiting regimes pre- 
viously described, the results are also valid for a free 
quantum particle on a torus if one uses the mapping 
{L — ► t — if; j — ► S-} (where t and if represent the fi- 
nal and initial time coordinates of the particle) . 

Our calculations demonstrate the role played by the 
geometry (curvature and shape) of the interface in con- 
trolling the size of the chain. 

As the size of the confining surface decreases, the poly- 
mer size is determined by the parameters of the confin- 
ing surface instead of chain length. The crossover be- 
tween the free Gaussian chain limit and the confined limit 
is controlled by a characteristic geometrical localization 
area A. A is proportional to a 2 for the cases of spheres 
and cylinders, and there are two characteristic areas pro- 
portional to a 2 and b 2 for the case of curved torus. For 
the cone, the localization area is variable and is deter- 
mined by the position of the singular point (the vertex) 
about the chain and by the angle a. Explicit closed for- 
mulas are derived for the dependence of the mean-square 
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FIG. 1. The geometrical data that characterize the sphere (A), the cylinder (B), the cone (C) and the torus (D): a is the 
radius of the sphere, of the cylinder and of the circular cross-section of the torus; 6 is the radius of the axial circle of the torus; 
a is the angle of the cone at the vertex. 

FIG. 2. The mean-square end-to-end distance ((R — R') 2 ) normalized to LI for a Gaussian polymer chain embedded on a 
cone, as a function of t — and a. 

FIG. 3. The mean-square end-to-end distance ((R — R') 2 ) normalized to LI for a Gaussian polymer chain embedded on a 
cone, as a function of t — and a: 3-D representation. 

FIG. 4. The mean-square end-to-end distance ((R — R') 2 ) normalized to a 2 for a Gaussian polymer chain embedded on 
different surfaces, as a function of j£ and t — The symbols joined by thin lines represent the data for the polymer on the 
surface of the curved torus in R 3 : thick solid line, the cylinder in R 3 ; o — o torus, t = 10; □ — □ torus, t — 8; O — O torus, t = 6; 
A— A torus, t = 4; V— V torus, t = 2; thick dotted line, the sphere S D_1 in R D . 
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Fig. 1 (Brownian Motion and Polymer Statistics on Certain Curved Manifolds 
Mondescu & Muthukumar) 



2 (Brownian Motion and Polymer Statistics on Certain Curved Manifolds - Mondescu/Muthukumar) 
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Fig. 4 Brownian Motion and Polymer Statistics on Certain Curved Manifolds: Mondescu/Muthukumar 
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